home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / ddisna.z / ddisna
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDDDDDIIIISSSSNNNNAAAA((((3333FFFF))))                                                          DDDDDDDDIIIISSSSNNNNAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DDISNA - compute the reciprocal condition numbers for the eigenvectors of
  10.      a real symmetric or complex Hermitian matrix or for the left or right
  11.      singular vectors of a general m-by-n matrix
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DDISNA( JOB, M, N, D, SEP, INFO )
  15.  
  16.          CHARACTER      JOB
  17.  
  18.          INTEGER        INFO, M, N
  19.  
  20.          DOUBLE         PRECISION D( * ), SEP( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      DDISNA computes the reciprocal condition numbers for the eigenvectors of
  24.      a real symmetric or complex Hermitian matrix or for the left or right
  25.      singular vectors of a general m-by-n matrix. The reciprocal condition
  26.      number is the 'gap' between the corresponding eigenvalue or singular
  27.      value and the nearest other one.
  28.  
  29.      The bound on the error, measured by angle in radians, in the I-th
  30.      computed vector is given by
  31.  
  32.             DLAMCH( 'E' ) * ( ANORM / SEP( I ) )
  33.  
  34.      where ANORM = 2-norm(A) = max( abs( D(j) ) ).  SEP(I) is not allowed to
  35.      be smaller than DLAMCH( 'E' )*ANORM in order to limit the size of the
  36.      error bound.
  37.  
  38.      DDISNA may also be used to compute error bounds for eigenvectors of the
  39.      generalized symmetric definite eigenproblem.
  40.  
  41.  
  42. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  43.      JOB     (input) CHARACTER*1
  44.              Specifies for which problem the reciprocal condition numbers
  45.              should be computed:
  46.              = 'E':  the eigenvectors of a symmetric/Hermitian matrix;
  47.              = 'L':  the left singular vectors of a general matrix;
  48.              = 'R':  the right singular vectors of a general matrix.
  49.  
  50.      M       (input) INTEGER
  51.              The number of rows of the matrix. M >= 0.
  52.  
  53.      N       (input) INTEGER
  54.              If JOB = 'L' or 'R', the number of columns of the matrix, in
  55.              which case N >= 0. Ignored if JOB = 'E'.
  56.  
  57.      D       (input) DOUBLE PRECISION array, dimension (M) if JOB = 'E'
  58.              dimension (min(M,N)) if JOB = 'L' or 'R' The eigenvalues (if JOB
  59.              = 'E') or singular values (if JOB = order. If singular values,
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDDDDDIIIISSSSNNNNAAAA((((3333FFFF))))                                                          DDDDDDDDIIIISSSSNNNNAAAA((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              they must be non-negative.
  75.  
  76.      SEP     (output) DOUBLE PRECISION array, dimension (M) if JOB = 'E'
  77.              dimension (min(M,N)) if JOB = 'L' or 'R' The reciprocal condition
  78.              numbers of the vectors.
  79.  
  80.      INFO    (output) INTEGER
  81.              = 0:  successful exit.
  82.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.